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I. Introduction 

In the asymmetric simple exclusion process (ASEP) on the integers Z a particle 
waits exponential time, then moves to the right with probability p if that site is 
unoccupied (or else stays put) or to the left with probability q — 1 — p if that site is 
unoccupied (or else stays put). A formula was derived in [9] for Py(xg(t) < x), the 
probability distribution function for xe(t), the position of the £th particle from the 
left at time t, given the initial configuration Y = {yi < y2 < • • •}■ The set Y may be 
finite or infinite on the right. The formula will be stated below. It is a sum over all 
(finite) subsets S of Y of integrals of order \S\. 

For step initial condition, when Y = Z + (the positive integers), three pleasant 
things happen: 

(i) For each k one can sum over all S with |,S| — k, thus replacing a sum over all 
finite S C Z + by a sum over k. 

(ii) There is a combinatorial identity that replaces the integrands by simpler ones 
which are symmetric in the integration variables. 

(iii) The resulting A> dimensional integrals are coefficients in the expansion of a 
Fredholm determinant. 

This led to a representation of the probability distribution in terms of a Fredholm 
determinant [10J, which was amenable to asymptotic analysis [TTj W\ |8| IT]. 

With step Bernoulli initial condition Y is not deterministic, but each site in Z + 
is, independently of the others, initially occupied with probability p. In this case also 
we have (i)-(iii), and this made asymptotic analysis possible [T2l IT]. 
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In another direction, Lee [5] considered the cases when Y = mZ + , and found that 
(i) and (ii) hold when m = 2, but only (i) holds when m > 2. Even for m = 2, 
(iii) does not seem to hold since the fc-dimensional integrals are not, or not obviously, 
related to a determinant or PfamanQ 

Nevertheless it seems worthwhile to consider a generalization, which we call 
periodic step Bernoulli initial condition. Here we have an m-periodic function n — > p n 
from Z to [0, 1] and assume that initially a site n G Z + is occupied with probability 
p n . The most that can be expected in this generality is that (i) holds, and we shall 
show that it does. 

Here are the formulas alluded to and the result of the present paper. We assume 
q 7^ 0, set r = p/q, and recall that the r-binomial coefficient [^] is defined as 



(l-r)(l-r 2 )---(l-r £ ) 

We define 

and then ^ 

i( x , k, o = i{x, k, . . . , &) = J] /(6, ^) n 

All indices in the products run from 1 to k. Notice that I(x,k,£) depends on t, 
although it is not displayed in the notations. 

Finally, given two sets of integers U and V we define 

a(U, V) = #{(«, v) : u eU, v e V, and u > v.}. 

Theorem 5.2 of [9], slightly modified, is the formula 

k 

scy JCr JC-r i = \ 

\S\>£ 

where k = \S\ and 

S = {sx,...,S k }, (Si < S 2 < • • • < Sfc) 
Qifc = (-l)^^-l)/2 r ^-D/2-M 



Earlier work on TASEP with initial condition 2Z+ or 2Z is in d [3] . 
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The sum is taken over all (finite) subsets S of Y with \S\ > I. The contour Cr is the 
circle with center and radius R, which is assumed so large that the denominators 
P + <z£i£j ~~ £i are nonzero on and outside the contours. (All integrals are given the 
factor l/2m.) 

When Y is not deterministic the probability P(xe{t) < x) is obtained by averaging 
the right side of ([[]) over all Y . The focus of the present note is (i), the computation 
of the sum over all S C Y with \S\ = k, and then the average over all initial configu- 
rations Y C Z + . What will be important here are only those ingredients of (CQ) that 
depend on Y and S. These combine as 

i=l 



For step initial condition the sum of (j2J) over all S equals 0, p. 838] 



r 



k 

fc(fc+l)/2 JT |_ 



= i 6 • • • Cfe - t 



fc-z+1 ' 



For step Bernoulli initial condition the sum over S C Y followed by the average over 
initial configurations Y is [T2l p. 830] 

k 

fc(fc+l)/2 TT P /q\ 

For y = mZ + the sum of (jSJ) over S 1 equals [5] 

*(fc+l)/2 TT * (A) 



To state the formula we derive here for periodic step Bernoulli initial condition 
we define 

<p(i,n) = 1 - p n +PnT k ~ l+1 

and m x m matrices A{ with entries 



.4 



1 



{iv<k) m iV 'P* 11 <P{i,n) iip<v, 

n=fi+l 



{^■■■ik) m - EI <P& n) 

n=l 



X < 



v+m— 1 

n 



if // > v. 



(5) 
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The indices p and v run from to m — 1. 

The result is that the average over Y of the sum of (J2j) over S CY with \S\ = k 
is equal to r k<yk+l ^ 2 times 



;i o • • • o) JJ a 

V 1 / 



(6) 



This is, in words, the sum of the entries in the top row of the matrix product. In the 
product, matrices with lower index % are on the left. To recapitulate: 

Theorem. For periodic step Bernoulli initial condition the probability ¥(xe(t) < x) 
is given by the right side of (JTJ with J2scy re pl ace d by J2k>e anc ^ ^ ne P r °duct (j2J) 
replaced by r k( - k+1 ^ 2 times (JHJ). 

In the simplest special cases, p n = p when n = v (mod m) and p n = otherwise. 
If we take < ^ < m, then (jHJ) is equal to 

k 



i=2 



which we find is equal to 

A- 



n 



j ,*(&'"&)" , -(i-p + /"-<=-' + i; 

when v — and 

it 

(6 • • • tk) m - u ] [ 



f = ite---a.) m -(i-p + p^ m ) 

otherwise. Alternatively, it is given by just the last formula if instead we take 
< v < m. 

In the course of the proof we obtain an analogous statement for what we may 
call general step Bernoulli initial condition in which the function p : Z + — > [0, 1] is 
arbitrary. We define B{ to be the semi-infinite matrix with entries 

v-l 

& v p» n (pit,™) tip<v, 

Bi,»,v={ " =M+1 (7) 

if p. > v, 
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the indices staisfying < //, v < oo. We shall see that the average over Y of the sum 
of (j2D over S CY with \S\ = k is then equal to T k ( k + 1 )/ 2 times 



10 •••) 1 [ B, 



i=l 



fl\ 
1 



Observe that this is what is obtained by formally letting m — > oo in the periodic 
case. Observe also that if p takes only the values or 1 we are back in the case 
of deterministic initial configuration, and what we obtain can be seen to be just a 
restatement of (pQ). 



II. Proof of the theorem 



We assume first that p : Z + — > [0, 1] is arbitrary and want to compute the average 
over Y of the sum of (EJ) over S C Y with \S\ = k. As in [12] we do this in the 
opposite order. We first take a fixed S and average over all Y D S. (Later we will 
take the sum over S.) Also, as in [12] , we may assume that Y C [1, N], and at the 
end we let N — > oo. (In fact all we shall use is that iV > Sfc.) 

1. The average over Y 

The part of (j2J) that depends on Y is r a( - s ' Y \ while the probability of an initial 
configuration Y is 

Y[pn JJ (1-Pn), 

n£Y ne[l,N]\Y 

so we have to compute 

e ^n*. n (!-^)- 

ScYc[l,N] neY ne[l,N]\Y 

Any element of Y which is larger than Sk does not affect a(S, Y), so this sum may 
be written product of sums, 

v 5cyc[i, Sfc ] nev ne[i,s fe ]\y 

The second factor equals 1, since it is 

II (Pn + (l"Pn)), 

ne(s fe ,AT] 



e no. 

yc(s fc ,iv] «ey 



n < 

ne(s k ,N]\Y 



(9) 
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when written as a sum of products. 

It remains to consider the first factor in ([9]). We use the notation 

Yt = (si_i, n Y, (s = 0), 

so Y is the disjoint union S U Y\ U • - • U Y^. Since the number of elements of S greater 
than or equal to Sj is k — i + 1, we have 

k k 

a(S, Y) = a(S, 3)+^2<t{S, YJ = k(k + l)/2 + ^(k - i + 1) \Yi\. 

i=l i=l 

Therefore the first factor in (Q equals 

A- 



r^+^n^nf e ^ i,|y<i iiA. n ^-pS) 

neS i=l ^y<C(«i-i,s«) n£>i ne(s i _i,s l )\y i 

= n pn n ( e n n ^ - *o) • 



n e5 i=i v yjcCsi-i.si) nev» nG(s i _i,s l )\y i 
The sum over all Yj C (si_i, s») is 

n n (^ +i Pn+(i-p„)), 

ne(si_i,Si) nG(si_i,Si) 

as is seen by writing this as a sum of products. Hence ([9]) equals r k ( k+1 ^ 2 times 

k 

n n ( io ) 

nGS i=l n£(si_i,Si) 

The sum over S: general p 

This gives the average of T tT ( s,Y * over all Y D S. Now we are to multiply this 
by Yli £7% the second factor in (j2J), and sum over all S with [ jS' | = k. Since 5 = 
{si, . . . , Sfc}, the sum over all S equals 

e n U s *• n ( n ) 

0<si <•»<«*. i=l ^ ne(si_i,Si) ' 

In terms of the entries of the matrices given by (El) this equals 

oo 

si,...,s/ s =0 

which equals ©, the sum of the entries of the top row of the matrix product. (We 
used here that B iai)V = unless \i < v.) 
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3. The sum over S: periodic p 

We shall rewrite the sum ffTTj) using the periodicity of p. There are unique repre- 
sentations 

Si = Uiin + Vi, where Ui > and < Vi < m. 
We have Sj_i < Sj if and only if the following hold: 

if v i-i < Vi then u^i < Ui, 

if Vi-i > Vi then Ui-i < U\. 
If we set Ui = ti + • — h tj, so that s« = (ti + • • • + ti) m + t> i5 then the above becomes 

U > if t>i_i < Uj 

(12) 

^ > if Vi-i > Vi. 



With this notation the factor ^ Si in (fTTj) becomes 

/•-(*iH Hi) m—Vi 

the factor p Si becomes p Vi by the periodicity of p, and by the periodicity of (p(i, n) 
in n the last factor becomes 



ti m+Vi — l 



n *>< 



n=«i_i +1 



n n^M) ifv<_i<vi, 

jjj+m— 1 y to \ *i— 1 

n n^M) nvi-.i>vi. 



If we fix the Uj and first take the sum in ffTTT) over the t{ according to (TP21 we 
obtain 



1 



Vi — l 

(Ci • • • &) m Pm 11 <P{h n ) if < v i 

n=tij_i+l 



n=l 



X < 



Wi+m— 1 
n=tij_i+l 



if > fj. 



Then the product of this over i is to be summed over all t>, satisfying < v , < m. 
In terms of the entries of the matrices Ai given by (j^J) the sum equals 

TO— 1 

^ ^ -^l, 0,"»i -^2, di , t>2 ' ' -^-k, Vk—itVki 
vi,...,v k =0 



which equals ([6]). This completes the proof of the theorem. 
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